LETTER TO THE EDITOR 



Collective Excitations in Realistic Quantum Wires 



Arne Brataasf A G Mal'shukov| Vidar Gudmundsson§ and K A Chaof 

t Department of Physics, Norwegian University of Science and Technology, N-7034 
Trondheim, Norway. 

flnstitute of Spectroscopy, Russian Academy of Sciences, 142092 Troitsk, Moscow 
Region, Russia. 

§ Science Institute, University of Iceland, Dunhaga 3, IS-107 Reykjavik, Iceland. 



Abstract. We have used the Hartree-Fock Random Phase Approximation (HF- 
RPA) to study the interacting electron gas in a quantum wire. The spectra of 
intersubband spin-flip excitations reveal a considerable red shift with respect to 
single-particle HF energies. That signals on appearance of collective intersubband 
spind-density excitations due to the exchange interaction. The long wavelength 
dispersions of the intrasubband collective spin-density excitations are linear, but 
the sound velocities are renormalised due to the exchange interaction and screening. 
The in-phase intrasubband charge-density excitation has the long wavelength form 
q[— ln(q)] 1 ' 2 . We found good qualitative agreement of our results with experimental 
observations. 



A semiconductor quantum wire can be fabricated by applying a voltage with a 
microstructured gate to a 2D electron gas. The single-particle energy (SPE) spectrum 
typically consists of subbands separated by several meV. At a ID electron density 
about 10 6 cm _1 more than one subband can be occupied. In recent years progress 
has been made in spectroscopic study of such systems. In angular resolved Raman 
spectra of GaAs quantum wires |, §, §, collective spin-density excitations (SDE) and 
charge-density excitations (CDE, or plasmons) were observed. The measured spectra 
cover from low to high frequency and for low-frequency intrasubband excitations, the 
wave- vector dependence of the spin- wave energy was found to be linear. The correct 
interpretation of these experiments allows us not only to understand the interesting 
physical processes, but also to access important physical parameters. 

The exchange interaction is crucial to the collective intrasubband and intersubband 
SDE. Similar to the direct long range Coulomb interaction which leads to the 
depolarisation shift of single-particle excitations and to the appearance of collective 
plasma modes, exchange interaction gives rise to the red-shift of SPE. If the red- 
shift is sufficiently large, the collective SDE with a sizable oscillator strength splits 
off the continuum of SPE. For semiconductor quantum wells, such split-off appears 
in the Hartree-Fock Random Phase Approximation (HF-RPA) |, |, 0. Thus, it 
is important to perform a HF-RPA analysis on collective electron excitations in a 
realistic GaAs quantum wire with full Coulomb interaction, and compare the results 
with measured spectra[|[ [|. The selfconsistent and conserving^ HF-RPA is suitable 
for this task, because we will calculate the two-particle spectra but not single-particle 
properties. As in using any approximation, the HF-RPA calculation also contains 
error. However, our HF-RPA results agree very well with experimental measurements. 

In this Letter we will first outline the HF-RPA method to express the spin or charge 
correlation functions in terms of the corresponding spin- or charge-density induced 
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matrix, which satisfies an eigenvalue matrix equation. In a realistic quantum wire, 
for high energy intrasubband excitations and for excitation spectra when more than 
one subband is occupied, the self-consistent equations have to be solved numerically. 
However, for the case that only one subband is occupied, the analytical expressions of 
SDE and CDE are derived. With additional approximation applied to the HF-RPA 
analytical solutions, they reduce to known results. This may be coincidental, but is 
outside the central theme of this Letter. 

To define the quantum wire, let us start from an electron gas in a narrow quantum 
well with interfaces parallel to the x-y plane, and the well width shorter than any 
other relevant length scales. We consider the realistic experimental situation that 
only the lowest subband in the quantum well is occupied, so the motion of electrons 
in the z direction can be ignored. By applying a properly designed gate potential, a 
quantum wire along the x axis is fabricated with a longitudinal constant electrostatic 
potential, but a transverse parabolic one Vc(y) = m*Lo 2 y 2 /2 along the y axis, where 
to* is the electronic effective mass. The quantum wire has a finite length L x , and we 
assume periodic boundary conditions. In the absence of electron-electron interaction, 
the single particle eigenfunctions are simply ip n k(x,y)=Lx 1 ^ 2 e lkx ')p n k(y), with the 
corresponding eigenenergies E n k—huJo(n+l /2)+(h 2 k 2 ) / (2m*) , where n is the subband 
index, k— integer x 2n/L x , and ipnk(y) the harmonic oscillator wave function. The 
transverse confinement length of the electrons is lo=[h / '(to* 'ljo)] 1 / 2 . 

When the electron-electron interaction is turned on, we will use the HF 
approximation 

/ h * w2 , vf , , 2e 2 [.j^f \Mr')\\ , ^ 

{ 2^ V + Vc{y) + —J dr \ fb yzp[ }Mr) 

/ dr 2,h— TZ — T— ^a(r') = e a i>a{r) (1) 

to derive a complete orthonormal basis {i/j a (r)} of quasi-particle Hartree-Fock states, 
where e^ F is the Hartree-Fock energy, /„ is the Fermi occupation factor and k the 
dielectric constant of the surrounding medium. In terms of this basis set, we will 
use the corresponding time dependent HF (HF-RPA) to calculate the charge-density 
and the spin-density correlation functions, which describe the self-consistent linear 
response of the electron gas to an external perturbation. In this Letter, we restrict 
ourselves to the nonmagnetic ground state. 

For inelastic light scattering not close to the band gap resonance, the Raman 
intensities in polarised and depolarised scattering geometries are proportional, 
respectively, to the imaginary parts of the charge-charge correlation function and 
spin-spin correlation function^]. If we define p(q,t) as the Fourier transform of 
the charge-density operator and a(q,t) the Fourier transform of the spin-density 
operator along the spin quantisation axis, then, these correlation functions are 
X"{q, w) = -| / °° dte iuJt < | [i>(q,t),v(-q*, 0)] | > with v=a,p, where (| . . . |) denotes 
a thermodynamic average. In HF-RPA, they can be expressed as x iy (<f, u>) = 
^cd^cdi'l^ w )( c l exp — iq- r\d) in terms of the z/-density induced matrix K". The 
charge-density induced matrix satisfies 
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- - Vbc-MWcd&u)] . (2) 

cd 

where Va(, ;c d= J dr J dr'^*(T^Tpd(r)4 , b( r ')' l l ) c( r ')/\f h ~' r '\ ^ s the Coulomb matrix element 
in the basis of the HF single-particle states. When u=a, there is no direct Coulomb 
interaction (Hartree term) between the excited spins. If the exchange term (Fock 
term) is neglected, the spin-density excitation spectra are simply given by the quasi- 
particle Hartree-Fock energies. In HF-RPA, the spin-density excitations are shifted 
from the quasi-particle Hartree-Fock energies due to the exchange terms. As a result, 
collective spin-density excitation may appear in the spectra. This shift is a measure 
of the effective strength of the exchange term in the system. 

We will use the above HF-RPA analysis to investigate collective excitations in 
a multisubband system. This has to be done numerically, and the first step is to 
derive the quasiparticle basis set {ipa{r)} by solving (1) self-consistently via iteration. 
The functional basis must be sufficient large in order to ensure the required accuracy 
for the quasiparticle energy spectra and the spatial electron density of the ground 
state. When using this self-consistent quasiparticle basis set to calculate collective 
excitations, the required accuracy for the peak position of an excitation energy is the 
experimental linewidth ~0.1 meV. 

The intrinsic materials parameters of a GaAs quantum wire required in our 
calculation are to*=0.067too and k=12.4, where mo is the electron mass. The extrinsic 
sample parameters are the total electron density n, the transverse potential V c (y), 
and the length of the quantum wire L x . Our calculations are based on a sample 
with n=10.4xl0 5 cm -1 , which has been investigated experimentally in Ref. ||. The 
transverse potential V c (y) is determined by the experimental conditions. Since the 
experimental sample does not have a parabolic potential as can be seen from the 
intersubband CDE, we will adjust huj to fit the subband electron occupations. When 
we set hujQ— 7.9 meV (7o=120 A), the derived self-consistent Hartree-Fock subband 
spacing between the two lowest subbands is 5.4 meV at the zone center and 5.6 meV at 
the Fermi level. The electron densities in the two occupied subbands are no=6.4xl0 5 
cmT 1 for the lowest subband (labeled by i=0), and m=4.0xI0 5 cm -1 for the second 
lowest subband (labeled by i=l). The physical properties of the quantum wire are 
insensitive to its length, provided that the wire is sufficiently long. Therefore, we 
set L x = 1.0 fim. The experiment H was done at a temperature 1.7 K, which is the 
temperature used in our calculation. 

Let ki and Vi be, respectively, the Fermi wave-vector and the Fermi velocity of 
the ith subband. In the limit of long wavelength q—>0, the transverse excitation 
is forbidden (or allowed) if the excitation energy is low (or high), and so the 
corresponding physical processes are dominated by intrasubband (or intersubband) 
transitions. When q increases to fco-fci, the low energy intersubband transitions are 
activated. However, collective excitations with large energy and/or large q are strongly 
damped. To illustrate the main features of long wavelength excitations, let us consider 
a very simply case by neglecting the electron self-energies and the vertex correction. 
In this case we obtain analytical results for both the SDE and CDE energies. The SDE 
dispersion is linear in q, with sound velocities vq and v\ for intrasubband excitations 
in the two lowest subbands. With two subbands occupied, the CDE dispersion has 
two branches. The in-phase mode is 

Luj(q) = \qW2(v Q + v 1 )V(q)/h7r, (3) 

where V(q) is the Fourier transform of the Coulomb potential (e 2 / K)[(x-x') 2 +(y- 



4 



y') 2 ] 1 / 2 , and the out-of-phase mode is 

Up (q) = \q\\/v vi . (4) 

Such two-band SDE and CDE dispersions are exactly the same results as SchulzjlCj 
obtained for a two-subband Tomonaga-Luttinger model (TLM) . 

When we turn on the intrasubband and intersubband exchange interaction, as 
well as the screening, the SDE and CDE energies will be modified from the above 
expressions, and have to be derived numerically. In the polarised spectrum, the higher 
energy region is dominated by strong intersubband CDE. At q=0, the calculated 
CDE spectrum consists of a single peak which is characteristic to transverse parabolic 
confinement potentials [[[2). In the region of finite q the intersubband HF SPE energies 
show up at 5.5 meV. In the spin-flip depolarised Raman spectrum, we obtained the red 
shift of intersubband SDE with respect to HF SPE. The SDE has a dominant weight 
of the spectrum and appears at a resonance energy of 2.1 meV, which is shifted from 
the HF SPE, indicating the importance of vertex correction. These overall features of 
our numerical results have been observed experimentally^, [|. 

In Fig. [j] we show the calculated intrasubband SDE and CDE in the region of 
energy and frequency for which accurate experimental data are available. The mode 
with strong intensity is labeled as SDEO, and the next mode with weaker intensity is 
labeled as SDE1. The inset gives their dispersions, which are linear for small q, with 
the corresponding sound velocities 1^0=0. 9uq and iv,i=0.7«i. These sound velocities 
are smaller than the respective Fermi velocities because of the exchange screening of 
SDE induced by the intersubband virtual transitions, and the renormalisation due to 
the intra- and intersubband exchange interaction. The experiments^], ||| seem to have 
detected only u CT) o, which corresponds to the high energy SDE near the Fermi wave 
vector k of the lowest subband. It is probably the low electron density in the second 
lowest suband that causes the velocity of the SDE around k\ to be too small 
to be observed experimentally. The spin velocity v a fi agrees quantitatively with the 
experiment^]. For sufficiently large q, because of the finite length of our quantum 
wire, the Landau damping shows up in Fig. [I] in the form of enhanced intensities of 
satellite single-particle peaks around the SDE, instead of a broad band if the quantum 
wire is infinitely long. The mode SDE1 is not damped, since it does not enter the 
region of SPEs. 

The dispersions of the CDE and SDE are also plotted in the inset in Fig. 0. In 
the long wavelength limit, the in-phase CDE+ can be well fitted with q[— ln(<7)]~ 1//2 
according to as expected for ID electron gas(l(], The out-of-phase CDE- 

mode, corresponding to (|]) with an expected linear dispersion appears in our 
calculation as a weak band with a sound velocity 1.2(vo^i) 1 / 2 . This happens to be the 
same as the sound velocity of the SDEO, an accidental result for this specific sample. 
The CDE- mode has a much lower intensity compared to the CDE+ mode. Checking 
against the experiment^], we believe the the lower energy band observed in the low 
frequency polarised Raman spectrum is our calculated CDE- mode. The decay of the 
plasmon (CDE+) and the CDE- at higher q due to the Landau damping is stronger 
than the decay of SDEO and SDE1, as shown in Fig. [|. 

After the above complete study of a realistic quantum wire, perhaps it is 
worthwhile to mention some surprising findings for the special case that only the 
lowest subband occupied at zero temperature, namely, the pure ID system including 
spin degrees of freedom. In this case we can drop the band index. For SDE, the 
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eigenvalue equation has the form 

Kl{q,u) = — fk :f~ fk+ f F ~ I dk'V{k-k')K^{q,Lo), (5) 

where V(k) = (2e 2 / K)Ko(qlo) for the quasi-lD Coulomb potential (e 2 / k)[(x- 
.t') 2 +/q] -1 / 2 , which has a logarithmic singularity (— 2e 2 / k) lnql as q<— »0. Since 
the Hartree self-energy is canceled by the potential energy due to the positive 
background charges, the difference in exchange self-energy at the Fermi energy, 
A x (q)=X x (k +q/2)-Z x (k -q/2) is A x (q) = ± dk> [V(k') - V(k' - 2k Q )}. We 
define y(k) by K%=(f k _ q / 2 -f k+q / 2 )y(h), and at fc=fc rewrite @ as 

q/2 

dk' [-V(k')y(k + k') + V(k' - 2k )y{-k + k')] 

-q/2 

= 2tt [Tiuj - hqv - A x (q)] y(k ) . (6) 

In the limit of long wavelength, by substituting y(ko + fe')^y(fco) and y{— kg + 
k')~y(—ko) into the above equation, we see that the contributions of Coulomb 
potential in the exchange energy and the vertex corrections, both behave as q\n(q), 
exactly cancelling each other. Hence, the q ln(g) behaviour of low-energy single-particle 
excitation due to the exchange self-energy is removed from the SDE and the plasmon 
energies, and so only the exchange terms at momentum 2fco are important. An 
equation similar to (6) for y k at k = — fco can also be derived. In the long wavelength 
limit, we find the SDE energy 



Ua(q) = \q\v oy /l-2g , (7) 

where go=V(2k$)/2hirvo. The sound velocity is reduced with respect to the bare Fermi 
velocity. We have also solved for the CDE and found the plasmon dispersion 



u p {q) = \q\v y/2V(g)/hn + 1 - 2g . (8) 

For a long-range Coulomb potential V(q)~\nq, the plasmon energy is not affected 
by the exchange interaction within the HF-RPA. These results (Q) and (^), derived 
with full Coulomb interaction, differ from those obtained by Schulz[|l(| for the TLM 
including a nonsingular backscattering matrix element at 2ko, which mixes the right- 
and the left-traveling modes. Nevertheless, in (||), if we replace the Hartree- Fock 
energy by the bare single-particle energy, and omit the singular part around V(q=0) 
in the vertex correction, then our results exactly reduce to Schulz's results for the 
SDE and the plasmon energies. Whether this finding, as well as (|J) and (||), is simply 
a coincidence remains to be clarified. 

To close this Letter, we should mention that according to our analysis, the 
frequency of the lowest energy SDE decreases as q increases towards fci, and then 
vanishes at a value of q close to 2ki. This behaviour suggests a low temperature 
intrinsic instability of the electron gas against the formation of spin-density waves 
(SDW), namely, the Peierls instability. When this instability emerges from our mean- 
field analysis, we must take into account the SDW long-range order in the ground 
state. On the other hand, the results of the TLM with backscattering predict the 
absence of SDW long-range order, but instead the appearance of a slowly decaying 
Wigner crystal at 4fc F @. We guess that the TLM with backscattering describes the 
real ground state better than HF-RPA. Hence, we believe that our analysis based on 
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a ground state without SDW long range order is a better approximation than that 
including SDW in the ground state. 
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Figure captions 



Figure 1. Intensity of intrasubband SDE and CDE as a function of the frequency 
shift for different wave vectors q. The upper part shows the dispersion of intrasubband 
CDE and intrasubband SDE. The quantum wire is specified in the text. 



